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ABSTRACT 

We present a new way of describing the flares from Sgr A* with a self-consistent 
calculation of the particle distribution. All relevant radiative processes are taken 
into account in the evolution of the electron distribution and resulting spectrum. 
We present spectral modelling for new X-ray flares observed by NuSTAR, together 
with older observations in different wavelengths, and discuss the changes in plasma 
parameters to produce a flare. 

We show that under certain conditions, the real particle distribution can differ 
significantly from standard distributions assumed in most studies. 

We conclude that the flares are likely generated by magnetized plasma consis¬ 
tent with our understanding of the accretion flow. Including non-thermal acceleration, 
injection, escape, and cooling losses produces a spectrum with a break between the 
infrared and the X-ray, allowing a better simultaneous description of the different 
wavelengths. We favour the non-thermal synchrotron interpretation, assuming the in¬ 
frared flare spectrum used is representative. 

We also consider the effects on Sgr A*s quiescent spectrum in the case of a density 
increase due to the G2 encounter with Sgr A*. 

Key words: Galaxy: center - Galaxy: nucleus - accretion discs - black hole physics 

- MHD - radiation mechanisms: thermal - relativistic processes - methods: numerical 

- galaxies: jets. 


1 INTRODUCTION 


1.1 Sagittarius A* 


Sagittarius A* (Sgr A*) is the name given to the bright radio 
source of our Galactic Center. It was discovered in 1974 by 
|Balick fe Brown (1974 \ using the Green Bank 35 km radio 
link interferometer of the National Radio Astronomy Obser¬ 
vatory. Stellar motion around the non-thermal radio source 
shows that Sgr A* is highly compact (smaller than 0.01 pc 
i.e. 3 x 10 11 km) and th at th e stars orbit around a poin t mas s 
of 4.3 ± 0.5 x 10 6 Mq ( Eisenhauer et aT~|| 2005 Melia|[2007 1 . 
The stellar orbits provide the strongest evidence yet for a su- 
permassive black hole located in the center of our Galaxy at 
a distance of 8.3 ± 0.35 kpc ( Reid|1993 Schodel et al.|2002 


Ghez et al. 2008 Gillessen et al.| 2009i. Super-massive black 

holes (SMBHs) of millions to billions of solar masses are be¬ 
lieved to exist in the centre of most galaxies. Sgr A*, in our 


own galaxy, is the closest and best studied SMBH, making 
it the perfect source to test our understanding of galactic 
nuclei systems in general. But among all galactic nuclei that 
we have observed so far, Sgr A* has the peculiarity of being 
very faint in all wavelengths. Even though it may have been 
more active in the past (|Revnivtsev et al. 2004; Zubovas fe| 


Nayakshlnj|2012 Ponti et al. ||2012 1, today Sgr A* is one of 

the most under-lumin ous SMBHs we know , it is very faint 
with Lboi — 10~ 9 TEdd I Narayan et al. 19981 and accreting at 


a very low rate. The accretion rate has been constrained by 


polarisation measurements, using Faraday rotation (Aitken 
|et al.|2000{|Bower et al.|2003||Marrone et al.|2007 ) and is es 
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timated to be in the range 2 x 10 _a < M < 2 x 10~'M@ yr" 
Theoretical work suggest that Sgr A* is most likely accreting 
at the lower range of this interval |Moscibrodzka et al.[2009 
Drappeau et al.[ 2013 ). Dibi et al. (20121 have shown that for 
accretion rate above or equal to 1 x 10 -8 solar masses per 
year, the cooling losses become important in the modelling 
of the accretion flow and the resulting spectrum. This re¬ 
sult means that Sgr A* is the only black hole source known 
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2 Salome Dibi et al. 


where cooling could still be treated separately as a first ap¬ 


proximation. Along the same lines, Yuan et al. (2004) had 
shown that flare events, as these observed in the Galactic 
Centre, could not be detected if the accretion rate increases 
by a factor 10 from its actual value, because synchrotron self 
Compton emission from thermal electrons would increases 
substantially. This would explain why Sgr A* is the only 
source known to exhibit flaring activity. Also, |Yusef-Zadeh | 
|et al.|2009| are reporting observational evidence for IR flaring 
activity inversely proportional to the flux density. 


1.2 Multiwavelength observations 

The faint emission from Sgr A* has been observed in dif¬ 
ferent wavelengths giving us a broad band spectrum of this 
object from the radio to the X-ray (see reviews by Melia 


|fc Falcke||2001| |Genzel et al.|[2010| |Morris et al.||2012[ and 
references therein). From a few GHz up to 100 GHz, the 
radio spectrum extends as a rough power law F v oc v a with 
0.25 ^ a ^ 0.33. Above 100 GHz there is evidence for a 
millimeter/sub-millimeter (sub-mm) excess over the power 
law, extending almost to 1000 GHz. The nature of this ex¬ 


cess was discussed by Serabyn et al. (19971 and Falcke et al. 


(1998) who excluded the possibility of dust emission. The 


size of photosphere is predicted to be the smallest around 
this wavelength of 1.3 mm, where the excess is observed. 


And the black hole horizon, or its shadow (Falcke et al. 2000 
|Dexter et al.|20lo| ) could be detected for the very first time 
in the near future with new, very long base interferometry 
facilities such as the “Event Horizon Telescope” (Doeleman 
|et al.|2008|[2009|[Flsh et al.|201lD . 

Important progress has been achieved in the mid- 


infrared (MIR) to near infrared (NIR; e.g., Genzel et al. 
|2003| |Ghez et al.||2005| |Schodel et al.|[2011[ ) and sub-mm 
domains. But in the optical and in the ultra-violet wave¬ 
lengths, the Galactic Center is heavily obscured by gas and 
dust with 30 magnitudes of visual extinction. This obscur¬ 
ing medium becomes partially transparent to the X-rays at 
energies above 2 keV. Indeed, Sgr A* has a quiescent X-ray 


luminosity of a few 10 33 erg s _1 (Baganoff et al. 2003) which 
is about 10 11 times lower than the Eddington luminosity. 

Sgr A* is quite variable and we observe fast activity 
(bursts or flares) in the infrared and X-ray band emissions. 
A few times a day, Sgr A* experiences rapid increases in the 


NIR flux (Hornstein et al. 2002 Genzel et al. 2003 Gliez 

et al. 20041 Eckart et al. 

2006, 2008; Yusef-Zadeh et al. 2008 

Dodds-Eden et al.|2011 

Haubois et al. 

2012 

), where brighter 

1 are often asso- 

flares (> 10 mjy; Dodds-Eden et al. 

2011 


ciated with simultaneous X-ray flares (e.g. Baganoff et al. 

to 

o 

o 

h- 1 

Goldwurm et al.||2003; Porquet et al. 2003; Belanger 

et al. 

2005) Porquet et al. 2008j Nowak et al. 2012; Neilsen 

et al. 

20131. The typical timescale for such events is few 


thousand seconds, suggesting a common localized origin of 
the flares. The radio and sub-mm emissions are more sta¬ 
ble, i.e. they show less variability than the X-ray and NIR 
(see for instance Marrone et al.|2008 for the sub-mm flares, 


and |Yusef-Zadeh et al.|2 0rC)|for a study on the IR - sub-mm 
anti-correlation). The flat radio emission is most likely the 
synchrotron emission originating from an outflow of Sgr A*, 
the lower the frequency, the further away we are in the out¬ 
flow. So the radio wavelength, as well as the quiescent X-ray 
emission are coming from extended regions around Sgr A*, 


while the sub-mm, MIR, and flaring X-ray emissions orig¬ 
inate from a region very close to the SMBH. This second 
region is the one we are interested in and we explore in this 
paper. 

The MIR and NIR emission has been observed by the 
VLT and Keck (e.g. [Dodds-Eden et al.||2011 Schodel et al. 


|2011[ |Bremer et al.||2011| |Witzel et al.||2012| ). The X-ray 
variability has been observed by XMM-Newton, the Chan¬ 
dra X-ray observatory (e.g. Baganoff et al. 200 1| ) and also 
by Swift ( |Degenaar et al.|2013| ). Many new X-ray flares have 
been detected recently thanks to the Chandra 2012 Sgr A* 
X-ray Visionary Project]^]. From this 3-Ms campaign, 39 X- 
ray flares are reported, lasting from a few hundred seconds 
to approximately 8 ks, and ranging in 2-10 keV luminos¬ 
ity from 10 34 ergs s' to 2 x 10 35 ergs s 1 ( Nowak et al.||2012 


Neilsen et al.|[20l3| ). The new telescope NuSTAR ( |Harrison| 
et al.|2013 1 has released recently new X-ray flares data that 


we are using in our study (Barriere et al., submitted). Those 
show that the 3-80 keV emission is compatible with a pure 
power-law spectrum. 


1.3 Flare models 

The fast variability indicates that the origin of the flares is 
as close as few gravitational radii from the SMBH. How¬ 
ever, the nature of the physical processes responsible for the 
flares is still an open question. Different mechanisms have 
been proposed such as events of magnetic reconnection or 
other acceleration processes (e.g., Markoff et al.|20'oT Yuan 
et al. 2003 Liu et al. 2004| 20061, infall of gas clumps or 
disruptions of small bodies (Cadez et al. 2006 Tagger & 
Melia; 2006 Z ubovas et al.|201 2), adiabatic expansion of hot 


plasma or hot spot models (Yusef-Zadeh et al. 2008; Brod¬ 
erick & Locb 2006). By modelling the physical conditions 


around Sgr A* and fitting the observational data, we also 
aim at giving a possible interpretation of the phenomenon. 

Several studies have been devoted to the modelling 
of SgrA* flares. Some models include a precise descrip¬ 
tion of the flow geometry. For instance, the emission from 
SgrA* was interpreted in the framework of radiatively in¬ 


efficient accretion flows (Yuan et al. 2003). In this model, 


the matter properties (density, temperature, etc.) are com¬ 
puted through hydrodynamical equations including radia¬ 
tive losses. These properties depend on the distance to the 
black hole and each ring contributes differently to the overall 
spectrum. The outer parts of the accretion flow contribute 
significantly to the X-ray luminosity in the quiescent state 
(Quatacrt 2002; Baganoff et al. 2003) with only about ten 
per cent of the quiescent X-ray flux coming from the cen¬ 
tral part we are modelling (Wang et al.|2013 Neilsen et al. 
20131. 


However, even in models where the geometry is dealt 
with accuracy, most of the emission originates from the very 
central parts of the accretion flow, both in the quiescent 
sub-mm and NIR bands, and in the flaring sub-mm to X- 
ray bands. Moreover, the typical time scale of a flare (few 
thousand seconds, depending on the flare) is of the order of 
the orbital period at the inner most stable orbit of SgrA*, 
pointing again to a flaring region of only a few gravitational 
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Exploring plasma evolution during Sgr A * flares 3 


radii (tg = GM/c 2 ). Therefore, most attempts to model 
the sub-mm to X-ray spectrum of SgrA* in the quiescent 
and flaring states (excluding the radio emission) implicitly 
assume that the emission originates from a single homoge¬ 
neous, isotropic zone characterized by only few parameters 
such as the average electron temperature and density, the 


magnetic field intensity (e.g. Dodds-Eden et al. 2010 Liu 
et al.|2006 1 . Here we use the same approach. 


Most models agree on the the fact that the small emit¬ 
ting region is weakly magnetized (< few hundred Gauss) 
and faint (< 10 8 particles/cm 3 ) which now appear as stan¬ 
dard values (Dibi et al.|2012| Moscibrodzka et al.|2009| Dex 


|ter et al.|2009| ). These values are supported by observations 

that constrain the accretion rate via Faraday rotation, to a 
maximum of M ~ 10~' solar masses per year. As a simple 
check, taking this higher accretion rate limit and a “typical” 
bulk velocity at one or two gravitational radii of 10 % of the 
speed of light (as simulated in GRMHD models of Sgr A*), 
then p m ax — M/{A-kR 2 v) ~ 10 8 particle/cm 3 . 

Whatever the details of the accretion flow and the ra¬ 
diative processes responsible for the emission, the emitted 
spectrum depends drastically on the particle distribution. 
For the sake of simplicity, all models so far have assumed 
pre-determined particle distributions (Maxwellian, power- 
law, broken power-law, or a combinations of them) which are 
described by few parameters. The precise shape of the par¬ 
ticle distributions depends on the radiative and acceleration 
processes and can deviate significantly from the assumed 
ones. The present work aims at dealing more precisely with 
particle distributions. Fitting arbitrary distributions to the 
data is not possible with current coverage and sensitivity 
of the instruments. Rather, the shape of the particle distri¬ 
bution can be computed self-consistently with a Boltzmann 
equation assuming a physics described by a few parameters. 
Such an approach is common in the modelling of the high 
energy emission from X-ray binaries and other AGN (e.g. 
McConnell et al7||2002| |Rogers et al.||2006| |Belmont et al.| 


20081 but has not been applied to SgrA* yet. In this paper, 


we present spectra obtained by solving simultaneously an 
equation for particles and an equation for photons to pro¬ 
duce self-consistent particle distributions and spectra. These 
spectra are compared to broad band data of SgrA* to put 
constraints on the flare properties. 

This paper is organised as follows: In Section 2 we 
present the microphysics and numerical method. In Section 
3 we present the results and solutions for the quiescent and 
flaring spectra from Sgr A*. We study the plasma behaviour 
in two kinds of models, namely in systems where matter is 
trapped in the emission region, and in systems where matter 
flows in and out of the emission region. In Section 4 we end 
with our conclusions and outlook. 


2 METHOD 

The goal of this study is to model the plasma around Sgr A*, 
and in particular the particle distributions and the resulting 
spectra. In the following, we will note v the frequency of 
photons, 7 the particles Lorentz factor , and p = — l ) 1 / 2 

the particle momentum. 

We use the belm code (Belmont et al.! 2008|. This nu¬ 
merical tool solves simultaneously coupled kinetic equations 


for leptons and photons in a magnetized, uniform, isotropic 
medium of typical size R. In all models presented here, 
this size is set to R = 2 cg = 1.3 x 10 12 cm based on 
the size derived from the flare time scale variability (where 
rG = GM/c 2 is the gravitational radius). 

The implemented microphysics includes radiation 
processes as self-absorbed radiation, Compton scatter¬ 
ing, self-absorbed bremsstrahlung radiation, pair produc¬ 
tion/annihilation, coulomb collisions, and prescriptions for 
particle heating/acceleration. 


2.1 Radiative processes 


Synchrotron radiation is produced by charged particles spi¬ 
raling around magnetic field lines. It depends on the mag¬ 
netic field B whose intensity is described by the magnetic 
compactness: 


Ib = 


<ttR B 

m e c 2 8-7T 


(1) 


where m e is the electron mass, c is the speed of light, or is 
the Thomson cross section, and R is the size of the emission 
region. Synchrotron emission at frequency v from a single 
electron with momentum p 
sivity js{p,v) in erg s _1 Hz 


characterised by the 

emis- 

Ghisellini et al. 

1988 

1998 


Katarzynski et al. 2006 1 . Synchrotron emission typically pro¬ 


duces soft photons and cools the high energy emitting par¬ 
ticles. The cooling time of relativistic particles emitting at 
frequency v is: 

t synch = 1.29 x 10 12 x u~ 1/2 x B~ 3/2 ( s ) 

For typical values of B, we have 

t = 1.3 (iVl0 18 Hz)“ 1 / 2 (B/100G )“ 3/2 s 


( 2 ) 


This corresponds to very short time scales, and only parti¬ 
cles emitting at frequency lower than 10 12 Hz cool on time 
scales comparable or longer than the duration of a typical 
flare (1000s). These particles are not observed to contribute 
much to the total emission. Moreover in our study we are not 
modeling the emission below 10 12 Hz that is extended radio 
emission from the outflow. Low energy particles can also 
absorb photons through the synchrotron process. Such ab¬ 
sorption is described by the absorption cross section cr s (p, v) 
(Crusius & Schlickeiser 1986 Ghisellini et al. 1998| |. The 
joint effect of high energy particle cooling and low energy 
particle heating tends to thermalize the particle distribu¬ 


tions. It is called the synchrotron boiler effect (Ghisellini 


et al. 19881. 


Photons of the emission region can also be scattered by 
Compton interactions. The scattering of isotropic photons of 
energy hvo by isotropic particles of energy Eo = 70 m e c 2 is 
characterized by the resulting distribution of scattered pho¬ 
tons o c (po,uo — > u). The belm code uses the exact, Klein- 
Nishina cross section (Jones 1968 Belmont 2009). In the 


case of SgrA*, soft photons are up-scattered by high energy 
particles which brings them to high energy. This also cools 


the scattering particles. From equation (37) of Piran (20041, 
the typical inverse Compton cooling time is: 

x B- 1 ^ (*) 


tc 


= 3.1 x 10 10 x „- 1 / 4 


(3) 


Again, this time scale is much shorter than the flare dura¬ 
tion. 
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4 Salome Dibi et al. 


The effect of self-absorbed bremsstrahlung radiation is 
also computed. However for the inner most regions of the 
accretion flow, bremsstrahlung emission is a negligible com¬ 
ponent of the resulting spectra and it will not be discussed 
on this paper. 

Photon-photon pair production and pair annihilation 
are also implemented in the code. However, we aim at mod¬ 
elling the emission from SgrA* only below 100 keV where 
these processes are negligible. They were disabled in order 
to reduce the computation time. 

Rather than computing the path of photons out of the 
emitting region with Monte Carlo simulations, photons pro¬ 
duced in-situ are assumed to escape with a probability rep¬ 
resentative of the geometry. This probability depends on 
the photon energy. For instance, high energy photons do 
not inverse Compton scatter and can escape freely when 
the optical depth is large, low energy photons can be scat¬ 
tered so much that they remain trapped in the system much 
longer. We use the escape rate from |Lightman fe Zdziarskij 
(19871; Coppi (2000) that reproduces well the results of 


Monte Carlo simulations in a spherical geometry. At low en¬ 
ergy, synchrotron and bremsstrahlung can absorb photons 
before they escape. This modifies the escape probability in 
this energy range. We include the corresponding modifica¬ 
tions to escape probability derived from Sobolev ( 1957|) (see 
also Poutanen fe Vurm|2009L 


2.2 Particle acceleration and heating 


The particle distribution depends on the above mentioned 
radiative processes and on several additional processes. 

Coulomb collisions tend to thermalize the particle dis¬ 
tributions. The belm code include Coulomb cross sections 


derived from Nayakshin & Melia (1998). However, the very 


low density inferred for SgrA* make this process very inef¬ 
ficient. In all the results shown in this paper, real Coulomb 
collisions are completely negligible. 

In order to account for the observed high energy radia¬ 
tion, particles needs to be heated/accelerated to high energy. 
Solving for the particle distribution thus requires to address 
also the physics of particle acceleration/heating. Many pro¬ 
cesses have been proposed to account for high energy parti¬ 
cles (viscosity, reconnection, shocks, first and second order 
Fermi processes, etc.). However, the precise process at work 
in SgrA* is still unknown. Moreover the physics of these 
processes is not constrained well enough to have a precise 
modelling for their effect on the particle distribution. Only 
stochastic acceleration by MHD waves can be implemented 
easily in a Boltzmann equation for the particle distribution 
(see Li u et al.|2006 for an application to SgrA*). However, if 
particle escape is slow, it forms a quasi-Maxwellian distribu¬ 
tion and does not reproduce hard non-thermal distributions 
such as the one observed by NuSTAR. If particle escape is 
efficient, it can produce power-laws only if the accelerating 
rate has the same energy dependence as the radiative cool¬ 
ing, which is very unlikely ( Katarzyriski et al.|2006 l. There¬ 
fore we use very general, ad-hoc prescriptions, inspired from 
what is done for the corona of accreting black holes (such 


as Eqpair, Coppi 2000 ; or belm). We use two different 
channels to provide energy to the particles. 


• We mimic thermal processes by computing Coulomb 


collisions with a virtual population of hot protons (with 
temperature ksT p = 40 MeV). Real collisions are very inef¬ 
ficient and do no provide significant heating, whatever the 
proton temperature. Rather, this prescription aims at repro¬ 
ducing the effect of anomalous processes (such as viscosity) 
on the lepton distribution. Therefore, the heating efficiency 
is renormalised by an arbitrary constant so that to inject 
power Lth (erg s _1 ) into the emitting region. In the follow¬ 
ing, this free parameter will be described by the compactness 
parameter l t h = orTth/ (Rm e c 3 ). Such prescription not only 
heats the global distribution of particles. It also thermalises 
it. As the efficiency of the virtual collisions is enhanced, the 
efficiency of the associated thermalisation is also enhanced 
to an anomalous level. Anomalous heating is a common fea¬ 
ture of accreting systems, so such heating is not surprising 
even though the origin is debatable. 

• We model non-thermal processes by constantly inject¬ 
ing particles with a power-law distribution N( 7 ) oc 7 _s . 
This distribution is characterised by 4 parameters: the slope 
s, the minimal and maximal energies 7 m i n and 7 max respec¬ 
tively, and the normalisation. As we want this process to 
keep the number of particles constant, the re-injected par¬ 
ticles are taken from the lepton population itself, with a 
uniform probability, independent of their energy. In the fol¬ 
lowing, the minimal energy of the power-law will be set to 
7 min = 50 , so that particles are accelerated from the bulk of 
the distribution. Indeed, the thermal peak of the spectrum 
(around 10 12 Hz) imply an electron temperature around 
10 11 K. And the maximal energy of accelerated particles 
will be set to 7 max = 4.6 x 10°, large enough to reproduce 
the NuSTAR data. The high energy cutoff has not been 
confirmed by NuSTAR observations (Barriere et al., sub¬ 
mitted), and the possible physical processes responsible for 
the non-thermal component (turbulent acceleration, recon¬ 
nection, weak shocks) can accelerate electrons to very high 
energies. The normalisation is computed so that the non- 
thermal process injects into the region a power L nt h, de¬ 
scribed by the free parameter l nt h = <JTL n th/(Rm e c 3 ). The 
slope is also a free parameter of the model. 

Such prescriptions compete with all other processes to pro¬ 
duce complex distributions of particles. 


2.3 Modelling the particle dynamics 

Although observational evidence clearly indicates that the 
sub-mm to X-ray emission originates from a very small re¬ 
gion close to the black hole, the dynamics of the particles is 
very uncertain. Free, relativistic particles can travel though 
the emitting region in a light crossing time: 

t = R/c = 43 s (4) 

which is much shorter than the flare duration. 

However, the medium is magnetised so that particles are 
not free to move on straight trajectories. Instead, they are 
bound to the magnetic field lines. For magnetic intensity of 
1-100 G, even X-ray emitting particles have gyro-radii orders 
of magnitude smaller than the emitting region. Therefore, if 
the medium is turbulent and the magnetic field tangled, even 
the highest energy particles can be considered as trapped in 
the main flow. 

The detailed structure of the accretion flow and in par- 
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ticular the accretion velocity are not known. Therefore we 
investigate two extreme scenarios. 

2.3.1 The closed system approximation 


and the slope of the power-law s, and the injection compact¬ 
ness Zinj ■ The particle density is no longer a free parameter 
and results from the balance between injection and escape. 
When injected particles have a relativistic temperature, the 
steady-state optical depth is: 


On the one hand we consider that the accretion velocity is 
very small. In that case, particles remain a very long period 
of time in the emitting region. Radiative and acceleration 
processes set steady distributions of particles and spectra 
before particles escape from the system. 

This model is characterised by 5 free parameters: the 
lepton density n e related to the Thomson optical depth r 
by r = n e OT-R, the magnetic field compactness Ib defined 
by Eq. □ the power of the thermal heating and non-thermal 
acceleration characterized by the compactness parameters 
Ith and Z n th respectively, and the slope of the non-thermal 
heating process s. In this model without particle escape, the 
particle distribution results from the balance between ther¬ 
mal heating, non-thermal acceleration, and radiative cool¬ 
ing. At high energy, thermal heating is inefficient, and the 
particle distribution results from the balance between ra¬ 
diative cooling and non-thermal acceleration. In good ap¬ 
proximation, Compton and synchrotron processes have the 
simple cooling laws shown in Eq. [3] and [2] As acceleration 
tends to produce an electron power-law distribution of in¬ 
dex s, the steady distribution is also a power-law with index 
s' = s + 1. When synchrotron radiation is the dominant pro¬ 
cess, this produces a synchrotron spectrum of spectral index 
a = s'/ 2. At lower energy, the physics and the shape of the 
particle distribution are more complex. 

We solve this model numerically for different parameter 
sets presented in Figure |T| [3] and [3] 


2.3.2 The open configuration 

On the other hand, we also consider the extreme case where 
matter flows in and out of the emitting region with an ac¬ 
cretion velocity approaching the speed of light. In that case, 
particles escape the emitting region on time scales compara¬ 
ble to the light crossing time, i.e. comparable also to the ra¬ 
diative times scales. Escape can therefore compete efficiently 
with radiation and acceleration processes. This model will 
be referred to as the open configuration. The distribution of 
the matter entering the emitting region needs to be given 
Amj (y). We assume that non-thermal acceleration occurs 
only in the emitting region, so that particle entering this 
region have a thermal distribution described only by 2 pa¬ 
rameters: its temperature 0i n j = fcslinj/(m e c 2 ) and its nor¬ 
malisation. The former is set to 8i n j = 13 in all models. The 
latter is described by the injection compactness: 



where h\ n j is the total number of particles injected into the 
emitted region per unit time, and the last equality holds for 
thermal distributions with relativistic temperature (#j n j >> 
1). Once in the emitting region, particles are assumed to 
escape on a typical time scale t es c = R/c, which corresponds 
to an escape probability p e3C = R/{c t esc ) = 1. 

This model is described by 4 free parameters: the mag¬ 
netic field compactness Ib , the non thermal compactness Z n th 


1 ^inj 
47T #inj 


( 6 ) 


In that model, the steady particle distribution results 
from the balance between injection and non-thermal acceler¬ 
ation versus both escape and radiative cooling. In this con¬ 
figuration, the shape of the steady state distribution is more 
complex than in the closed system. At high energy, particles 
cool before they escape. As in the closed system, the lep¬ 
tons form a power-law distribution of index s' = s + 1, and 
emit a power-law synchrotron spectrum of spectral index 
a = s'/ 2. At low energy, particles escape before they cool 
and the steady state electron distribution is a power-law 
of slope s' = s. This produces a synchrotron spectrum of 
spectral index a = s/2. The particle distribution and pho¬ 
ton spectrum thus exhibit a break, whose energy depends 
on the relative efficiency the cooling and escape. As far as 
synchrotron radiation is the dominant cooling process, the 
break in the photon spectrum is: 


Vbreak = 2.97 X 10 


s — 1 
3.6- 1 

R 


3 x 10 12 


III 

R/c , 


(7) 


Such a model was for instance proposed by |Dodds-Eden| 
et al. (2010) to explain the flaring X-ray emission without 


violating the NIR upper-limits. They used a broken power- 
law distribution. However, depending on the parameters, 
processes other than synchrotron emission can contribute 
to the physics. Also, a self consistent cooling break is not 
sharp and extend over a significant frequency range. Here 
we extend their conclusion by solving self-consistently for 
the particle distribution. 


3 SGR A* RESULTING SPECTRA 

The flare duration varies but the typical time is about 3000 
seconds. The radiative time scales and the thermalisation 
time are much smaller than the flare duration. The particle 
and photon distributions are in quasi steady state at each 
moment of the flare, and the flare evolution is directly gov¬ 
erned by the evolution of model parameters, here namely 
the acceleration processes described here by the parameters 
lnth and 1th- Therefore we will mostly present and discuss 
the results from steady state simulations. I.e. we aim at re¬ 
producing separately the quiescent and flaring spectra by 
changing the value of only few parameters. Doing so, we 
can say that the flaring state is just a transition governed by 
the increase or decrease of few physical parameters. What is 
driving these changes is subject to interpretation, but know¬ 
ing what needs to be modified should give us some good first 
insight into the physical processes at work. 

Table 1 summarises the characteristics of our different 
models. 
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Table 1. Description of simulations 


Model 


Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


m 


spectrum 

configuration 

B 

n e 

Ith 

^nth 

^inj 

s 

Pesc 

tk/f-b 

T e 

Lbol 



Gauss 

part/cm 3 

XlO -5 

XlO -5 

XlO -4 




XlO 10 K 

X 10 36 erg s -1 

Quiescent 

Closed 

154.3 

4.6 x 10 6 

2 

1 

0 

3.60 

0 

0.15 

11.4 

1.4 

Flare 

Closed 

48.8 

4.6 x 10 6 

7 

3 

0 

2.40 

0 

6.32 

35.0 

4.8 

Flare 

Closed 

48.8 

4.6 x 10 6 

5 

4 

0 

2.20 

0 

5.79 

33.4 

4.3 

Quiescent 

Open 

175.3 

3.3 x 10 6 

0 

1 

4.64 

3.60 

1 

0.08 

7.84 

1.0 

Flare 

Open 

175.3 

3.3 x 10 6 

0 

9.8 

4.64 

2.28 

1 

0.09 

7.90 

3.8 

Flare 

Open 

175.3 

3.3 x 10 6 

0 

11.4 

4.64 

2.13 

1 

0.09 

7.87 

4.8 

Flare 

Open 

34.5 

1.4 x 10 s 

0 

100 

200 

2.60 

1 

97.16 

7.79 

5.8 

Quiescent 

Open 

175.3 

9.9 x 10 6 

0 

1 

14 

3.60 

1 

0.24 

7.86 

3.8 


The first column gives the Figure’s number. The second one gives the state of the spectrum we are trying to reproduce (in 
quiescence or during a flare) and the configuration of the model (closed or open region). B is the magnetic field which is a free 
input parameter. n e is the density which is also an input parameter in the closed configuration cases (Figure 00 and 0 but 
an output of the simulations in the other cases. is the prescription for thermal heating, and Z nt h the prescription for non- 
thermal heating/acceleration . Z; n j is the injection of particles for the open configuration cases, s is the spectral index used in the 
non-thermal prescription for acceleration. P eS c gives the probability for particles to escape from the system, e^/ei, is the ratio of 
kinetic energy over magnetic energy (8rr f mc 2 (7 — I) /B 1 ) and is computed after the different models have run. T e is the 

temperature of the thermal part of the spectrum. Because our lepton distributions are calculated self-consistently we never have 
a perfect thermal distribution, but this temperature corresponds to the closest match between the real distribution and the pure 
Maxwellian (plotted in dotted lines in all our distributions as a comparison), this is of course an output of our models as well. 
Finally, the last column gives the total luminosity of the spectrum. In the closed configuration case, this luminosity is the direct 
result of the thermal plus non-thermal compactness parameters, while in the open configuration it results from the balance between 
injection, escape, and cooling. We remind that in all our models we have R = 2rQ, 'tmin = 50, 7 max = 4.6 X 10 5 , and 0, nJ = 13. 


3.1 Data 


Sgr A*’s SED is made of different observations that are vari¬ 
able and in most cases have not been observed simultane¬ 
ously in different wavelengths. We have chosen a set of data 
representative of the overall variation of the spectrum. In all 


figures, the black radio points are from Falcke et al. (19981 
and Zhao (2003), the red radio points are recent ALMA ob¬ 


servations from Brinkerink, Falcke et al. submitted). The 
black far IR upper limits are from Serabyn et al.| (19971 


and Hornstein et al. (20021, the green MIR data are from 
Schodel et al. (|2011 1, the pink NIR lower point (in the qui¬ 


escent spectra) and the cyan NIR upper point (in the flare 
spectra) are from Ghez et al. (20041; Genzel et al. ( 2003), 
and |Dodds-Eden et al. (20111. The green “bowtie” is from 


Bremer et al. (2011) and is one of the few slopes that has 


been observed so far in the IR. Several flare observations 
seem to be consistent with this value of the NIR spectral 
index around — 0.6 ± 0.2 (Ghez et al.||2005 Gillessen et al. 


2006 Hornstein et al.|2007 Bremer et al.||20TTL 

The black “bowtie” in the quiescent X-ray is from 
|Baganoff et al.| ( |2001| |2003| ) and is an upper limit for the 
central emission because it is contaminated by thermal 
bremsstrahlung from the outer accreting matter. The or¬ 
ange “bowtie” is a Chandra flare from |Nowak et al. ( |2012[ ). 
Finally the blue data points (dark and light blue) are two 
flares observed with NuSTAR on July 21st and October 17tli 
2012 respectively ( [Barriere et al. 2014J. Even-though the X- 
ray flares may seem different in shape and slope, they are 
both acceptably fit by an absorbed power-law, and their 
photon indices are not significantly different (2.23 +q' 22 and 
2.041 q'2o for the July 21st and October 17th flares, respec¬ 
tively). Barriere et al. 2014, investigated the presence of a 
cutoff in the October 17th flare, but found that it is not re¬ 
quired by the data. The other wiggles in this spectrum (one 


may see a ” V” shape in the low energy part of the spectrum) 
are not significant either. One need to keep in mind that the 
error bars show the 1-sigma confidence range, which means 
that an acceptable fit does not need to go through them all 
but 3 out of 9. 

Three X-ray flares are shown on the first flaring spectra 
(Figure 0 and 0 , but then we consider only the modelling 
of one of the NuSTAR. flares that extend to higher energies. 
For a better reproduction of the spectrum, we need to move 
to the “open configuration” where we present some possible 
spectra for the July flare or the October flare. 


3.2 Sgr A* spectra from a closed region 

Figure [l] shows a spectrum for the quiescent state of Sgr 
A*, together with the emitting steady state lepton distri¬ 
bution. For this first spectrum, we consider a density of 

4.6 x 10 6 particles per cubic centimetre (which corresponds 
to r = 4 x 10 -6 ) and we keep this density to study the 
closed region, i.e. we consider that the number of particles 
is kept constant in the quiescent and flaring states. In qui¬ 
escence, the magnetic field is 154.3 Gauss, the plasma is 
magnetically dominated with tk/tb = 0.15. The thermal 
heating is twice the non-thermal one and corresponds to 

9.6 x 10 35 erg s -1 and 4.8 x 10 35 erg s _1 respectively, so 
that the total emission reaches 1.4 x 10 36 erg s _1 in qui¬ 
escence. We can see in the resulting spectrum in Figure [T] 
that the non-thermal component is not dominant, and this 
is not only due to the low value of l n th but also because 
of the steep injected slope s = 3.6. Nevertheless, this non- 
thermal component is important in order to reproduce the 
lower NIR data point. The thermal part contributes mainly 
to the sub-millimetre bump. In this case, we can notice how 
the thermal part of the lepton distribution differs from the 
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Figure 1. Quiescent spectrum from Sgr A* (top panel) and the 
associated lepton distribution (bottom panel) in a closed system 
configuration. On the spectrum, the black data points are taken 
from \Yuan et af .j \2003\ with the X-ray “bowtie” corresponding to 
an upper limit for the quiescent state of Sgr A *. The data points 
on the spectrum are described in details in section 3.1. The blue 
curve component of the spectrum corresponds to the synchrotron 
process, while the red one corresponds to the Compton process. 
The Bremsstrahlung contribution is not visible in the scale of this 
plot. On the electron distribution, the solid line is the shape of 
the calculated distribution from which the spectrum comes from, 
while the dotted lines indicate a pure Maxwellian plus power-law 
components for comparison. 

standard Maxwellian shape on the bottom panel of Figure 
[I] Indeed, for particle energies around p = 10 2 me, the differ¬ 
ence between the calculated distribution and the standard 
shape in dotted line, can reach almost two orders of magni¬ 
tude. The steady state particle distribution is sharper than 
a pure Maxwellian. Above 7 = 100, synchrotron cooling 
overcomes the anomalous thermalisation, and the distribu¬ 
tion cuts more sharply than a pure thermal one. This also 
produces a sharper sub-mm bump as is illustrated in the 
resulting spectrum (top panel of Figure [I]). On Figure [ 2 ] we 
plotted the spectral shape resulting from the dotted line of 
Figure [l] We can see that the quiescent spectrum would be 
much wider, reaching the far-IR upper limits. The novelty 



Figure 2. Quiescent spectrum from Sgr A* resulting from 
the “standard” distribution consisting of a simple Maxwellian plus 
power-law (dotted line in the bottom panel of Figure [7]j. 


of our work is illustrated by the difference between Figure [l] 
and[ 2 j that results from the careful and detailed calculation 
of the lepton distribution. 

Starting from similar conditions as in the quiescent 
state of Figure [I] Figure [3] shows a spectrum for the flaring 
state of Sgr A*, together with the lepton distribution. The 
spectrum is dominated by synchrotron self-Compton emis¬ 
sion (red line), even-though the non-thermal synchrotron 
(blue line) has a non negligible contribution to the total 
spectrum. The emitting region is the same as in the quies¬ 
cent state with the same density of particles. However the 
magnetic held has dropped from 154.3 to 48.8 Gauss, the 
plasma being now kinetically dominated with tkflb = 6.3. 
This dramatic change could be interpreted as being due to 
magnetic reconnection, a physical process that could be at 
the origin of the haring event. When the magnetic held is 
rearranged, some magnetic energy is converted to kinetic 
energy, thermal energy, and particle acceleration. In this 
way we have a decrease of the magnetic held strength and 
an increase of the two parameters 1 th and l n th represent¬ 
ing the thermal and non-thermal acceleration respectively. 
In this case, to have the Compton dominated spectrum in 
the Haring state (Figure |3|, the thermal heating increases 
from 9.6 x 10 35 erg s -1 to 3.4 x 10 36 erg s -1 , and the non- 
thermal one from 4.8 x 10 35 ergs -1 to 1.4xl0 36 ergs -1 . The 
non-thermal component has also a much hatter distribu¬ 
tion in the Haring state, meaning that the high energies are 
more populated, while s is steeper during quiescence. With 
this model, during the hare, the total luminosity reaches 
4.8 x 10 36 erg s -1 . 

Figure [4] shows another potential hare model, together 
with its lepton distribution. This spectrum is similar to Fig¬ 
ure [3] except that the non-thermal synchrotron component 
is more important than the Compton one. The radius of 
the emitting region is the same as well as the density of 
particles. The magnetic held magnitude is also 48.8 Gauss. 
The difference comes from the balance between thermal ver¬ 
sus non-thermal heating. For this non-thermal synchrotron 
dominated spectrum, the non-thermal contribution l n th is 
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Figure 3. Flare spectrum from Sgr A* (top panel) and the asso¬ 
ciated lepton distribution (bottom panel) in a closed system con¬ 
figuration. The data points on the spectrum are described in the 
beginning of section 3. The blue curve corresponds to the syn¬ 
chrotron process, while the red corresponds to the Compton pro¬ 
cesses. The Bremsstrahlung contribution is too small to be visible 
on this scale. On the lepton distribution, the full line corresponds 
to the actual calculated distribution, while the dotted line is a 
standard Maxwellian plus power law as a comparison. In all our 
models (except for Figure [Ify , the spectra result from the calcu¬ 
lated particle distribution (full line), while the theoretical fixed 
distribution (dotted line) is just shown as an illustration. 


Figure 4. Flare spectrum from Sgr A* (top panel) and the as¬ 
sociated lepton distribution (bottom panel) in a closed system 
configuration. The data points on the spectrum are described in 
the beginning of section 3. The blue curve corresponds to the 
synchrotron processes, while the red corresponds to the Compton 
processes. It is the same as Figure [~] but for the case of a dom¬ 
inant synchrotron component with respect to the Compton one. 
On the electron distribution, the full line is the shape of the ac¬ 
tual distribution from which the spectrum comes from, while the 
dotted lines are pure Maxwellian plus power-law components as 
a comparison. 


more important than previously with a value of 1.9 x 10 3e 
erg s -1 and a slope of 2.2. The total luminosity is similar to 
the previous case with L=4.3 x 10 36 erg s _1 . 

The quiescent state is very well reproduced by this 
closed region configuration model, the sub-millimetre bump 
is clearly fitted by synchrotron emission which extends to 
the lower part of the variable MIR and NIR data, and we 
are not violating the X-ray limit represented by the black 
“bowtie”. According to the new results from [Wang et al. 
(20131 and Neilsen et al. (20131 saying that the inner region 
is dominated by non-thermal emission from combined weak 
flares that can contributes to 10% of the quiescent X-ray, we 
could be even too low in the X-ray luminosity (lower than 


10% of the observed flux). The flaring spectra are somewhat 
more marginal because the sub-millimetre contribution is 
too high compared to two black data points that are upper 
limits and end up below the spectra. On the other hand, the 
sub-mm part of the spectrum is also variable on the order 
of 20% and considering that we don’t have perfect simul¬ 
taneous data, it still provide a close enough interpretation, 
meaning that we are still within 20% of the actual data point 
values. In the MIR, the flaring spectra shown in Figures [3] 
and [4] are in the right range of luminosity, and we can also 
reproduce the X-ray flare fluxes. The NuSTAR (blue) and 
Chandra (orange) flare slopes are respected, while the trend 
of one of the MIR flare (green “bowtie”) is not well repro¬ 
duced at all. We have to keep in mind that our data are 
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not simultaneous and slopes in the MIR have been observed 
only few times, but still in this case our slope seems to be 
in contradiction with this observation. 

Comparing models with observations, the “a” pre¬ 


scription from Shakura & Sunyaev (19731 is still used to 


parametrize turbulence and quantify the angular momen¬ 
tum loss mechanism, and the process whereby gravitational 
binding energy is converted into radiation. The best physical 
interpretation of this a parameter is given by the mechanism 
of magneto-rotational instability (MR.I). For instance, Haw¬ 


ley et al. (19951 have performed three-dimensional magneto¬ 


hydrodynamic numerical simulations of an accretion disk to 
study the nonlinear development of the MRI, they obtained 
that the time average of a is 0.6 for the vertical field runs. In 
the study of advection-dominated accretion and black hole 
event horizon, Narayan & McClintock (2008) argued that 


for an advection-dominated accretion flow, the theoretically 
expected value of a is 0.1-0.3. We have the estimated num¬ 
bers for a for the different cases studied here. Assuming 
Keplerian assumptions, which is obviously a very simplistic 
approximation close to the black hole but allows us to check 
roughly that the orders of magnitude are not inconsistent 
with the first order a approximation, the viscous heating Q 
is related to the “a” parameter by: 

n 3 / GM \ 1/2 

Q= 2 aP {lp) 


( 8 ) 


where P is the pressure. Using the dimensionless constants 
from our model, the viscous parameter is derived from the 
following formula: 


1 Itn r 1 ' 2 

2 r e e 


(9) 


with t=2tg, t = n e aTR, O e = kT e /m e c 2 . We find that 
the viscosity parameters a resulting from models |T] [3] and 
[4] are 0.18, 0.18 and 0.16 respectively. These values are in 
good agreement with the theoretical predictions described 
above which illustrates that an anomalous thermal heating 
is common on the context of accretion disks. 




Figure 5. Quiescent spectrum from Sgr A* (top panel) and the 
associated lepton distribution (bottom panel) with thermal injec¬ 
tion and escape (open configuration). The data points are the 
same as in the previous quiescent spectrum on Figure [1] 


3.3 Plasma with escape and thermal injection 

For comparison within the open configuration, we want to 
reproduce the quiescent spectrum with the assumption that 
particles flows in and out of the emitting region. Figure [5] 
shows such a spectrum, which is similar to the simulated 
spectrum in Figure [T] This spectrum is a realistic and ac¬ 
ceptable solution for the quiescent state of Sgr A*. The mag¬ 
netic field magnitude is 175 Gauss with a resulting plasma 
density of 3.3 x 10 6 cm -1 that is also magnetically dominated 
with tk/tb = 0.08. We note that in this case, the thermal 
component does not really differ from the pure Maxwellian 
distribution (see bottom panel of Figure [5|. Indeed, in this 
range, radiative cooling is negligible, so that the balance be¬ 
tween the thermal particle injection and the uniform particle 
escape produces a steady state distribution that is almost 
the pure Maxwellian. The total luminosity of this spectrum 
is L = 1.0 x 10 36 erg s _1 , also equivalent to the previous 
quiescent fit. From this quiescent spectrum we next investi¬ 
gate the changes necessary in order to move to the flaring 
spectrum. 

Figure [6] shows a spectrum for a flaring state of Sgr 


A*, together with the lepton distribution. The flaring spec¬ 
trum is dominated by non-thermal synchrotron and repro¬ 
duce the NuSTAR July flare as well as an IR flare with a 
slope closer to the usually observed one (flat to slightly ris¬ 
ing in the power spectrum). As expected, cooling breaks are 
observed in the lepton distribution and in the photon spec¬ 
trum. These are not sharp but span at least one order of 
magnitude in frequency. The emitting region is the same as 
in the quiescent state and the magnetic held stays the same 
as well. The amount of injected particles is the same as in 
quiescence leading to a constant density. The only change in 
order to move from the quiescent to the hare spectrum is on 
the non-thermal component: the heating parameter l nt h is 
increasing by almost one order of magnitude, and the slope 
becomes hatter (from 3.6 to 2.3 during the hare) meaning 
that we have more particles in the higher energy part of the 
electron distribution. So, we must have some physical pro¬ 
cesses that accelerates the particles more efficiently in the 
haring state and creates a harder non-thermal distribution. 
As a consequence the total luminosity increases, reaching 
3.8 x 10 36 erg s" 1 . 
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Figure 6. Flare spectrum from Sgr A* (top panel) and the as¬ 
sociated lepton distribution (bottom panel) in the open configu¬ 
ration. The data points are the same as in the previous spectrum 
on Figure [3] The electron distribution shows also in dotted line, 
the pure Maxwellian and power-law curves as a comparison. 


Figure 7. Flare spectrum from Sgr A* (top panel) and the asso¬ 
ciated lepton distribution (bottom panel) with thermal injection 
and escape. The data points are the same as in the previous flare 
spectra on Figures |~j [7] and [7j The calculated electron distribu¬ 
tion (full line) and the theoretical one (dotted line) as a compar¬ 
ison. 


We can do the same exercise to reproduce the October 
NuSTAR flare. This is shown in Figure [7] for our best case 
scenario that is really similar to the model on Figure[6]with 
non thermal synchrotron with a cooling break responsible 
for the flare emission. This is not surprising as we explained 
earlier in the Data section 3.1, that both flares are not sig¬ 
nificantly different and should be modelled with a power-law 
shape as a fit. As for the July flare, the trigger of the event 
is on the non-thermal component of the lepton distribution 
that increases by a bit more than an order of magnitude 
with a prescribed slope of 2.1 which is a bit flatter than for 
the July flare. Beside that, the size of the emitting region, 
the magnetic field, and the density are the same as for the 
other flare and the same as in quiescence. For the Octo¬ 
ber flare model we have a slightly higher non-thermal power 
with a slightly flatter prescription for the acceleration, the 
total luminosity of this flare spectrum is 4.8 x 10 36 erg s -1 
We investigated an alternative scenario where the X- 
ray flares would be produced by synchrotron self Compton 
(SSC) emission, however we found that models that account 


for pure SSC as an emission mechanism have physical pa¬ 
rameters that are hardly compatible with what we know 
of the central region density. Moreover it leads to a more 
complex scenario where the large scale magnetic field and 
the density of the medium need to be modified during the 
flare event. Nevertheless, inverse Compton emission could 
still be a non-negligible component of the overall spectrum, 
especially at high energies, assuming a weaker magnetic field 
during the flare and a higher density medium. Figure [8] gives 
an illustration of some “power-law” shape X-ray emission 
that would be a combination of synchrotron and SSC. In 
this case the density has increased from 3 x 10 6 to 1 x 10 8 
cm' 1 and the magnetic field has dropped from 175 to 35 
Gauss moving to a kinetic dominated flow with e^/tb = 97. 

We think that the best model for the flaring state of Sgr 
A* is the one produced by non-thermal synchrotron with a 
cooling break as seen on Figure [6] and [7] because the trends 
of the multi-wavelength data are reproduced and only very 
few parameters need to be adjusted in order to move from 
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Figure 8. Flare spectrum from Sgr A* (top panel) and the asso¬ 
ciated lepton distribution (bottom panel) with thermal injection 
and escape. The data points are the same as in the previous flare 
spectra on Figures 000 and 0 The calculated electron dis¬ 
tribution (full line) and the theoretical one (dotted line) as a 
comparison. 


the quiescent to the flaring state. This is especially true if 
we consider that the green “bowtie” is a typical IR slope. 
The non-thermal synchrotron emission is a simple and ele¬ 
gant solution of the flaring event observed by Chandra and 
NuSTAR because the overall state of the medium does not 
change dramatically (for instance the density and magnetic 
field is kept constant). The acceleration of the electrons lead¬ 
ing to the more important and flatter non-thermal lepton 
distribution is the only modification, and this could be trig¬ 
gered by some plasma instabilities that are not modelled in 
details here. 

Even-though magnetic reconnections could also be the 
initial trigger, it can happen on very small scales and does 
not necessary lead to a drop of the global magnetic field 
magnitude. A possible sudden increase of the density (as in 
model0l can be interpreted as an accretion rate fluctuation, 
however such fluctuations of more than an order of magni¬ 
tude are most likely not happening every day in the Galac- 



Figure 9. NuSTAR X-ray light curve of the flare event of the 
21st of July 2012. The black data points is the unabsorbed flux 
between 3 and 79 ke V. When the detection was not strong enough 
to be significant, we have only plotted upper limits of three sigma 
(arrows). See \Barriere et al . | \20!J,\ for the observed light-curve. 
The red curve represents the X-ray light curve from our quiescent 
spectrum model in Figure[5\ to the flare spectrum model in Figure 

m 


tic Center and would be difficult to interpret. The NuSTAR 
data being consistent with a power-law shape to higher en¬ 
ergies points also in favour of the synchrotron scenarios as 
in models [6] and 0 

The sub-millimeter part of the spectrum is really stable: 
comparing the quiescent state on Figure [5] with the sub¬ 
millimeter part of the spectrum on Figure0 we have ex¬ 
actly the same contribution around 10 12 Hz. This is mainly 
due to the fact that the magnetic field is kept constant and 
the injected population is also constant. This configuration 
gives us in return a constant density. We have to keep in 
mind that the data are not simultaneous, nevertheless it is 
an interesting exercise trying to model several wavelength 
observations in the same time. In the future simultaneous 
observations are going to be very important for this kind of 
multi-wavelength study. 

We then looked at the time evolution between Figure [5] 
and 0 to reproduce the X-ray light curve of the July NuS¬ 
TAR flare. Using our self-consistent calculations, we can 
also model the time-dependent particle evolution in order 
to reproduce the flare light-curves. This approach has also 
been considered by Dodds-Eden et al. (20101 in the one zone 
cell approximation but for a given power-law distribution. 
The cooling time-scales being very short compared to the 
flare duration, the time evolution is entirely governed by 
the physics of the acceleration processes that are not clearly 
defined. Figure 0 shows the reproduction of the X-ray light 
curve between 3 and 79 keV for the same parameter setting 
as for Figure 0 and 0 During the flare, the non-thermal pa¬ 
rameter l n th evolves linearly with time from the quiescent 
value 1CF 5 to a maximum value, such that the averaged 
value over the flare duration is 9.8 x 10 -5 as in our flare 
spectrum 0 It reaches a maximum value at the peak, and 
decreases back immediately with a linear dependence. The 
slope of the accelerated particles is set to 2.28 during the 
flare event as in our flare spectrum 0 and to 3.60 in qui- 


© 2014 RAS, MNRAS POQ.fTIfTi] 



























12 Salome Dibi et al. 


escence. We note that before the flare event, Sgr A* is not 
detected by the X-ray satellite NuSTAR because it is too 
faint, and embedded in the diffuse/unresolved emission, in 
this case we simply plotted 3 cr upper limits. 


4 CONCLUSION AND OUTLOOK 


We are able to reproduce the quiescent spectrum of Sgr A* 
in two different scenarios: considering that the accretion pro¬ 
cess is very slow and that the same particles remain a long 
period of time in the emitting region, and considering that 
the accretion process is very efficient, with and accretion 
velocity close to the speed of light; so that particles only 
remain in the emitting region on short time scales compara¬ 
ble to the radiative time scales. To model the flaring state 
however, we favour the second scenario, that allows a bet¬ 
ter interpretation of the sub-millimetre and infra-red part 
of the spectrum. The flaring state spectrum is best repro¬ 
duced by a plasma that has the same low magnetic field 
as in quiescent, and the same amount of injected particles. 
More efficient non-thermal heating processes are responsi¬ 
ble for the flaring event, and a flatter non-thermal distri¬ 
bution of electrons is present. Besides this change, all other 
parameters stay the same when moving from the quiescent 
to the flaring spectrum (Fig ure |5 | and |6| . Our conclusions 
are in good agreement with (Dod d~Eden et al.| |2010l who 
also favoured non-thermal synchrotron processes and a cool¬ 
ing break in order to explain the observed IR and X-ray 
flares. However, in our study we do not make the hypoth¬ 
esis of magnetic reconnection as an energy power for the 
flares, and our conclusions do not favour this particular hy¬ 
pothesis. As in our best case scenario (Figures [5] and [6] or 
[Tj), the magnetic field is not required to drop significantly. 
An important drop in the magnetic field amplitude has also 
important consequences on the sub-millimetre and thermal 
part of the spectrum that we also model here, other param¬ 
eters have then to be carefully adjusted in order to maintain 
the sub-mm shape in reasonable values, so we think other 
acceleration mechanisms are more likely to be happening. 
Reconnection mechanisms could also occur in very localised 
regions, and particles would diffuse away from the reconnec¬ 
tion sites and radiate in a field which has not reconnected, 
so we would not notice any significant global drop of the 
magnetic field amplitude. In our study, we end up with a 
plasma density of 3.3 x 10 6 particles per cubic centimetre, 
which is a reasonable value according to observations and 
theoretical work. But what would happen to the quiescent 
spectrum (Figure [l] or [5| if the density increases by a fac¬ 
tor three as expected to happen now when the cloud G2 is 
falling into the Galactic Center? As reported by |Gillessen| 
|eTaL |20l2| ), a dense gas cloud approximately three times 
the mass of Earth is falling into the accretion zone of Sgr 
A*, but nothing noticeable has been observed yet from Sgr 
A*. Figure [lO] represents such a prediction, it has exactly 
the same settings as the model described on Figure [5] but 
the density is three time higher (we have more particle in¬ 
jection). The model predicts a flux increase in the sub-mm 
bump (10 12 10 13 Hz), however the current emission is not 

well constrained in this band. If the source stays in quies¬ 
cence, we do not expect a particular increase in the IR, and 
we have some emission in the ultra-violet due to the first 




Figure 10. Quiescent spectrum from Sgr A* with the same con¬ 
ditions as in Figure [5| but assuming an increase of the density 
by a factor three. All the observed data (of quiescent and flaring 
state) has been kept on the figure. The associated lepton distri¬ 
bution is shown on the bottom panel. 


Compton component that is unfortunately not detectable. 
Even in the X-ray, if the increasing density by a factor three 
does not trigger a flare event, we do not expect a significant 
increase from the quiescent X-ray level. Overall it could well 
be that we are not detecting any striking changes. 
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